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K'^ (or 7r~) in the K~tt^ (or tt^tt^) centre-of-mass (cm.) frame; (ii) angular 
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distribution of e" in the e~e^ cm. frame; and (iii) the correlation between the 

meson and lepton planes. We also derive CP-violating observables obtained 

by combining the above decays with the conjugate processes B K^Tr~e~e'^ 

and B — > it' w'^ e'^ . 

PACS number (s): ll.30.Er, 13.20.He 
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Typeset using REVT^iX 



I. INTRODUCTION 



The purpose of this paper is to investigate the angular distribution of the decays B 
K^n^e^e'^ and B Tr^vr+e^e^, when the Ktt and vrvr systems are the decay products of 
K* and p respectively. The aim is to derive the detailed consequences of the effective short- 
distance Hamiltonian describing the four-fermion interaction b s{d)e~e^. This work may 
be regarded as an extension of previous studies of the exclusive processes B K*{p)e~e'^ 
IHJ^ that were limited to the kinematical variables si (the invariant mass of the lepton 
pair) and cos^^ (the angular distribution of /~ in the cm. system). The additional 
information we provide is the distribution in cos6p, where 6p is the angle of the (vr^) 
in the K~n~^ (7r~7r^) cm. frame, and the dependence on the angle (p between the e~e^ and 
K^Ti^ or TT^TT^ planes. This information is sensitive to the polarization state of the vector 
meson K* (p), and thus provides a new probe of the effective Hamiltonian. 

An important aspect of the calculation is the prediction of CP-violating observables 
that can be obtained by combining information from B and B decays. These observables 
probe a term in the Hamiltonian proportional to Im (V^bKr^/^fe^s) (i^ the case of B K*) 
and Im (\4f,l^*^/V^{,VJ^) (in the case of i? p). While the numerical estimates of these 
asymmetries in the standard model turn out to be small, the formalism we present is also 
applicable to more general Hamiltonians transcending the standard model. The dependence 
on the variable (p can be an especially useful probe of CP violation, as has been demonstrated 
in the analogous case of 7r+7r^e+e~ 

II. MATRIX ELEMENT 

We are concerned with the matrix element of the decay B{p) P{ki)P'{k2)l^{qi)l^{q2)-i 
where PP' = K'tC^ or 7r^7r+ and I = e,p. Introducing the linear combinations ||5|J^ 

k = ki + k2, K = ki-k2, q = qi + q2, Q = qi - q2, (2.1) 

the short- distance Hamiltonian for b —>■ fl^l~ (/ = s or d) is 
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- 2cf fia^,^{m^PR + mfPL)hlYl 



where Pl,r = (1 =F 75)/2 and 



cf = -0.315, Cio = -4.642, 



(2.2) 



(2.3) 



with 



Cg = cg + (3ci + C2)<^ 5'(mc, Si) + A„ 5f(mc, s;) - g{mu, si) 



4.224 + 0.361 



;i + Xu)g{mc, Si) - Xugirriu, Si) 



Xu 



ub uj 

V V* ' 



+ 



(2.4) 



(2.5) 



g 8 4 2 

g{mi, si) = -- \n{mi/mb) + 1^ + gVi ~ g(2 + vdVl'^-Vi 



1 / 1 + , 

m , — tiT 



Q{yi — 1)2 arctan ■ 



(2.6) 



1 - VI - ViJ J ' ' Vvt - 1 

where = and ?/i = Amf /si. The elhpses in the above represent numerically insignificant 
terms involving the Wilson coefficients C3,...,cq (see, e.g., Ref. |0]). The corresponding 
matrix element is 



cf - c,,){P{h)P\h)\h,PLh\B{p)) 



2c^' 



,efr 



■{P{h)P\k2)\fia^,q^ (mtPR + tti^Pl) b\B{p)) 



ClO -ClO 



. (2.7) 



Observe that the coefficient df contains both a weak phase (associated with the imaginary 
part of Xu) and a strong phase [connected with the imaginary part of ginii, s;)], thus opening 
the way to CP-violating asymmetries between B and B decays. 

The hadronic part of the matrix element, describing the transition B PP', can be 
written in terms of -B — > ^ form factors {V = K* or p) 
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{P{h)P'{h)\h,PLb\Bip)) 
with the convention eoi23 = +1, 



(2.8) 



= Dvik')!^ - te,,.pK^p^kPg+{q') T \k, [^7+(g')A + q'g^{q') 

± k,{{q . W) [g+{q') + \q'Kq')\ + \i[9+{q')^ + ^Vl^')] } + •■•} 



(2.9) 



where the elhpses denote terms proportional to q^ which may be dropped in the case of 
massless leptons and 



A = (M|-M2), W^ = K^-ik^, ^^ Ml Ml, 

k'^ 



(2.10) 



Throughout this paper, we neglect the lepton mass, and assume the above form factors 
to be real, the absorptive parts due to real cc and uu states having been included in the 
functions g{mc,si) and g{mu,si) appearing in the short-distance coefficient cf^ [Eq. (|2.4| )]. 
Furthermore, we will limit ourselves to neutral B mesons, i.e. B = B^ . 
The function Dyik"^) appearing in Eqs. (|2.^ ) and (|2.9| ) is defined via 

487r2 



2\\2 



\Dy{k')\ 



-5ik'-M^), 



where we have used a narrow-width approximation, and 

Ai/2(P,M2,M|.,) 



(2.11) 



(2.12) 



with the triangle function 



A(a, 6, c) = a'^ + + c- 2{ab + bc + ac). 



(2.13) 



We adopt the B — > K* and B ^ p form factors g, /, h, g±, a+ given by Melikhov and 
Nikitin (see Appendix The matrix element ( p.7|) can then be written compactly as 

-MsD = ^^^^^^//{(^W^'^^Al + K^VL + k^ZL)lYPLl + {L^ R)}, (2.14) 



where 



XL,R = Dvik^){icf T cio)^(g^) ^{rrib + mf)g+{q 



(2.15) 



yL,R = --Dy{e){{cf T cio)/(g') - -^{m, - mj)[g4q')A + gV(g')]}, (2-16) 



zl,r = {q ■ W)z'L^ji - ^vl. 



(2.17) 



with 



r 9r'^^ r 1 i ^ 

Dv{k'){{cf T Cro)a4q') + -^{m, - m^) [g^{q') + -q^q')] }• (2.18) 



The matrix element for the corresponding antiparticle channel can be obtained from ( |2.14D 
by means of CPT invariance. In fact, we have 



(2.19) 



for the conjugate processes B — > K'^n e and B ^ n vr+e e"*", where x,y,z are related 
to X, y, z by 



(2.20) 



Here 0w and S denote the weak and strong phases respectively that appear in the matrix 
element A^sd- 

III. DIFFERENTIAL DECAY RATE 



Squaring the matrix element ( p.l4|) and summing over spins, we obtain for the decay 



B PP'l-l+ lOl 



\MsBiB ^ PP'l-l+)\ 

X |2Re {ylzL + y^^Zij) [{k ■ q){q ■ K) - {k ■ Q){K ■ Q) - si{k ■ K) 
+ 2Re {xlzL - x*^Zr) \{k ■q){k- Q){q ■ K) + SiSp{K ■ Q) ~ {k ■ qf{K ■ Q) 



- si{k ■ K){k ■ Q) 

+ 2Re {xRy*j, - x^yl) [siK\k ■ Q) - {k ■ Q){q ■ Kf + {k ■ q){q ■ K){K • Q) 

-si{k-K){K-Q) 

+ {\XL? + kill') [-2{k- q){k ■ Q){q ■ K){K ■ Q) - siSp{K ■ Qf + {k ■ qf{K ■ Qf 
- siK\k ■ Qf + {k ■ Qf{q ■ Kf + 2si{k ■ K){k ■ Q){K ■ Q) 

+ {\yL? + \yR?) [ - SiK' + {q ■ Kf - {K . Qf 
+ {\zl\^ + \zr\^) [ - SiSp + (k ■ qf - {k ■ Qf 
+ 2e^,o.pk^K''q'^QP [{K ■ g)Im {xlvI + ^rVr) + Im {y*pZR - vIzl) 
- {k ■ Q)Im {xIzl + x*pZr)^ |, 

with sp = k"^. Introducing the shorthand notation 

1/2 1 



X 



{k.qf-sisp ^-X'/\Ml,Si,Sp), 



(3.1) 



(3.2) 



we find that (m; = 0) 



k-K^Ml-Ml,^ isp, 



(3.3a) 



k-q = -{Ml - si- Sp), 



(3.3b) 



k ■ Q — X cos 9i, 



(3.3c) 



q- K ^ pXcos9p + ^{k-q), 



(3.3d) 



q-W ^ pXcosOp, 



(3.3e) 



K ■ Q = ^{k ■ Q) + (3 {k ■ q) cos6i cos 6p — {sispf^"^ sin 6i sin 6p cos0 



(3.3f) 



e^^.a^k^K-q'^Q^ = - (s^Sp)^^/?^ sin sin sine/., 



(3.3g) 
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K' = 2{M'p + M'p,) -sp = (e' - Hsp, Q' = -si. 



(3.3h) 



To write the differential decay rate in compact form, we define the auxihary functions 



F^L,R = (3[X'z'+{k-q)yLA 



(3.4) 



F2L,R — 



(3.5) 



[3X{siSpf'^XL 



(3.6) 



so that 



216^8^3 \^tbVt}\ pXI{si,sp,9i,6p,(j))dsidspdcos9idcos6pd(j), (3.7) 



with 



I = Ii + I2 cos 201 + J3 sin^ 9i cos 20 + J4 sin 29i cos (p + I5 sin 6'^ cos 4> + Iq cos 6'/ 

+ Ij sin 6'; sin + /s sin 29i sin + /g sin^ 9i sin 20. (3.8) 



The functions Ji, . . . , Jg have the following form (see Appendix ^ for details): 
^3 



^ ( I FsL P + I F2L n sin^ 9p+\F,l\' cos' 9p 



+ {L^ R), 



(3.9a) 



IuFslI' + IFslH sin^ 9p - cos^ 



(3.9b) 



/3 = (|i^3Lr - \F2Ln sin' 9p + {L^ R) 



(3.9c) 



I4 = Re {F^lF;^) sin 2^p + {L ^ R), 



h = 2Re (Fi^Fg*^) sin20p _ (L ^ i?), 



/e = 4Re {F^lFD sin^ 9p - {L ^ R) 



(3.9d) 
(3.9e) 
(3.9f) 



/r = 2Im {F^lF^l) sm20p - {L ^ R), 



(3.9g) 



Is = Im {F,lF*l) sin 20p + (L ^ i?), 

Jg = -2Im (F^lF;^) sin^ 9p + {L^ R). 
The physical region of phase space is defined through 

^ ^ {Mb - v/sp)^ (^P + Mp,f ^sp^ M^, 



(3.9h) 
(3.9i) 

(3.10a) 

0<(/)^27r, -l^cos^p^l, -l^cos^z^l. (3.10b) 

Note that 6i is the angle of the /~ in the cm. frame; 6p is the angle of tt^ (or K") in 
the 7r^7r+ (or K^ti^) system; is the angle between p^- x and p;- x p;+ (in the case 
of the 7r~7r+ final state) and between p^- x p^+ and p;- x p;+ (in the case of the K'n'^ 
final state), pj denoting the 3-momentum vectors of the corresponding particles in the B 
rest frame. The z-axis is chosen along the total momentum vector of the PP' system in the 
B rest frame. This definition will also be retained in the case of the antiparticle reactions 
B 7r^7r+e^e+ and B K^-K^e^e^ (with p;^+ replacing p^^- in the latter case). 

The amplitudes Fh^r {i = 1, 2, 3) defined in Eqs. (|3.4| ) -(|3.6| ) are closely related to the 
transversity amplitudes Aq,A\\,A±_ sometimes used in connection with the angular distri- 
bution of the four-body final state arising from decays of the form B —>■ V1V2 |[n|. The 
differential decay rate in this alternative notation is again given by Eqs. (|3.7|)-(|0|), with 
the functions FiL r written in terms of Aq, Ah, A± as follows: 



where the normalization factor 
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FlL,R - 
F2L,R = 
F3L,R = 

2i%77Vf3 



OL,R 



N ' 

A\\L,R 
A±L,R 



1/2 



(3.11a) 
(3.11b) 
(3.11c) 

(3.12) 



has been chosen in such a way that F = \Aq\^ + l^iip + \Ax_\ 



TABLE I. Values of the functions J-'" [Eq. (|^], using p = 0.19 and r/ = 0.35 



121 











•-'2 


B ^ TT TT^e 


76.2 


76.8 


25.4 


-76.8 


B K-TT+e-e+ 


261.4 


278.8 


87.2 


-278.8 



IV. ANGULAR DISTRIBUTIONS 

We derive from the differential decay rate, Eq. ( p.7|) , one-dimensional angular distribu- 
tions of interest, namely dT/dcosOp, dT/dcosOi, and dV /dip. These distributions, as well as 
the observables calculated in Sec. |V], depend on different combinations of the short- distance 
coefficients Cg^, cio and the form factors f, h, g±, a+. 



A. Decay rate as a function of cos 6p 



Integrating d^T over the variables si, sp, cos 9i, and 0, we obtain 



rfcos^P 2i%7 \ tb tf\ yi 3 
with 

Ji = Jt sin^ Op + cos' Op, J-'' = I ^i'^PXdsi dsp, t=l,2, (4.2) 

where ll''^ stand for the coefficients of sin' Qp and cos' Qp in the expressions for Ji 2 given 
in Eqs. ( |B1|) and (B2) of the Appendix. The absence of a term odd in cos6p is connected 
with the fact that the PP' system is in a pure L = 1 state. As a consequence, the forward- 
backward (FB) asymmetry in the PP' system, defined as 



—dcosUp — / —dcos&p 

d cos 9p J_i d cos 9p 



^ ^-1'----'-^ (4 3) 

f dV ^ ^ dV ^ ' 

/ ^ — dcost^p+ / —dcosdp 

Jq dcosUp J_i dcostip 

vanishes. The predictions for J'-''^ are contained in Table |. 
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TABLE II. Estimate of the functions Ki according to Eq. ([4.5|). 





Ki 


K2 


K3 






Ke 


K7 


Ks 


K9 


B ^ TT vr^e 


152.8 


-17.3 


-14.1 








-61.8 








-0.05 


B K-Tr+e-e+ 


534.4 


-69.6 


-64.6 








-179.0 








-0.4 



B. Decay rate as a function of cos di 

Integration of d^T over si, sp, cos 9p, and yields 

= ^15^ \V,.V,}\\K, + cos 2^, + K, cos^O, (4-4) 



where 



~ y liPXdsi dspdcosOp. (4.5) 

Our results for the parameters Ki are tabulated in Table ||. Observe that K4, K^, Ki, and 
Kg vanish in our model. The FB asymmetry of /~ in the l^l^ cm. system is 

Kr/2 -0.19 for 7r"7r+, 

^FB = Y^^^ = (4.6) 
Ki-K2/^ 1-0.16 for K~'K+. 

This result reproduces the FB asymmetry of the lepton calculated in previous analyses of 
the exclusive channels B K*{p)e~'e^ 0)0] • 



C. Decay rate as a function of cf) 

Finally, the distribution in the angle between the lepton and meson planes, after 
integration over other variables, takes the simple form 



dV GlM%a^ 



{K, - ^K2) + ^(ifa cos 20 + Kg sin 20) 



(4.7) 



with the Ki^s shown in Table ||. 
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V. CP- VIOLATING OBSERVABLES 

We now focus on CP-violating observables that can be constructed by combining infor- 
mation on 13 and B decays, namely 

B ^ K-n+e-e+ and B ^ K+n-e-e+, (5.1a) 

or 

B — > 7:~^ e~ e'^ and B — > Tr^vr^e^e^. (5.1b) 



With the definition of cos6'p, cos9i, and given after Eq. (p.lO|) , the differential decay rate 
for B and B decays is 

d^T\^ = Afl{si, sp, Ou Op, (j))(3XdFS, (5.2a) 

d^T\B = ^fI{sl, sp, Oi, Op, (j))(3XdFS, (5.2b) 
where A/" is a normalization factor 

and dPS represents the phase-space element 

dPS = dsidspd cos Old cos Opdcf). (5.4) 

The function I is given in Eq. ( |3.8|) , whereas the function I is obtained from / by the 
substitution 

h,2,3Aji^L,R; yL,R] Zl,r) > h,2,3A,7{xL,R] yL,R, Zl,r) = A,2,3,4,7, (5.5a) 

hfi,8,9{xL,R] yL,R] Zl,r) > —hfi,%,%{^L,R-, yL,Ri ^L,r) = ~hfi,8,9, (5.5b) 

where x, y, z are defined via Eq. ( p.20[ ). 
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We now define the sum and the difference of the differential spectra d^T\§ and d^T\B as 
follows: 



Cardiff = ^^ 



(Ji - Ji) + {h - h) cos 261 + (Jg - h) sin^ 61 cos 20 + {h - h) sin 261 cos 
+ (/5 + Is) sin 9i cos + (Je + /e) cos 9i + (Jy - Jy) sin 9i sin + (Jg + Is) sin 29i sin 
+ (Jg + /g) sin^ sin 20] /3XdPS, (5.6) 



(Jl + Ji) + (/2 + COS 29i + (Jg + Jg) sin^ 9i COS 20 + {h + /4) sin 2^; COS . 



+ (Is - h) sin 9i cos + (Je - h) cos + (Jy + h) sin 6^; sin + (Jg - /s) sin 26'; sin ( 



+ (/g - Jg) sin^ 9i sin 20 



Recalling Eq. (|4.5| ) and defining 

1 



K, 



D 



Ml 



pXdsi dsp 



/•! 



lid cos 9p, 



-1 Jo 

we can then construct the following CP-violating observables: 



(^cp) = Y 
-'0 



(5.7) 



(5.8) 



(5.9a) 



{As) = Uk3-k,), 

-'0 



(5.9b) 



Jo 



(5.9c) 



-'0 



(5.9d) 



In 



(5.9e) 



(A) = ^(i^7''-^7''), 
-'0 



(5.9f) 



1 

7o 



(5.9g) 
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(5.9h) 



in which the quantity Iq is defined as 



Io = {K, + K,)--{K2 + K2). 



(5.10) 



Note that the asymmetries (^cp), (^3,4,7) involving the differences (Ji — li), i = 1, 2, 3, 4, 7, 
can be obtained from a measurement of the difference dTdm- On the other hand, the asym- 
metries (^5,6,8,9) require only a measurement of rfFsum and, in principle, can be determined 
even for an untagged equal mixture of B and B. 

The asymmetry (Aqp) is simply the asymmetry in the partial decay rates of -B — >■ 
K^TT^e^e'^ and B K^Ti^e^e^ (or B Tr^vr+e^e^ and B 7r^7r+e^e+). The other 
asymmetries represent CP-violating effects in the angular distribution of these processes. 
All of these effects have their origin in the CP-violating imaginary part of the coefficient 



In the case f = d, this imaginary part is of order rj, but in the case / = s it is reduced by 
an extra factor A^. In both cases, the weak phase present in the coefficient cf^ [Eq. ( p. 41) ] 
is further suppressed by a factor of order (3ci + C2)/cg ^ 0.085. All of these factors explain 
the very small magnitude of the asymmetries listed in Table |II[ The result for the partial 
width asymmetry (Acp) reproduces that obtained in previous literature 0. 

The small value of (Aq) implies that the FB asymmetry of the electron is opposite in 
sign for B and B decay [Eq. ( ^4.6| )j. This result agrees with the statement in Ref. [|1^. 



Our results for dT / d cos 6 p and dV/dcj) are new consequences of the short- distance 
Hamiltonian for h — » s{d)e^e^ , which test the polarization state of the vector meson in 
B — > K*{p)e~e~^ . The asymmetries (Ak) are CP-violating observables that can be obtained 




for f = d. 



(5.11) 



VI. CONCLUSIONS 
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TABLE III. Estimates of the average CP-violating asymmetries (Af^) in units of 10 ^ (10 ^) 
for the B K* (B p) transition. 

















(^5) 


(^6) 




(^8) 


(^9) 


B - 


^ K-TT- 

vs 




-e+ 


2.7 


-0.6 


-2.0 


5.2 


-4.6 





0.6 


-0.04 


B - 


K+TT- 




-e+ 


















B - 


VS 


'e~ 


-e+ 


-1.7 


0.1 


0.4 


-1.4 


1.2 





-0.1 


0.006 


B - 




'e~ 





















from a comparison of the angular distribution in B and B decays in the conjugate channels. 
Numerical estimates have been obtained within the framework of the standard model, but 
the formalism presented can be applied to other Hamiltonians. As seen from Table |ITT|, the 
asymmetries predicted for B K* are exceedingly small, and any significant effect would 
signal a non-standard source of CP violation. In the case of B p, asymmetries of 1-2% 
are predicted for (Aqp), (^5), and (^46). It should be noted that our predictions apply to K* 
(p) production on the mass shell. This implies that the PP' system is in a pure p-wave state, 



and explains why the distribution in cosOp given in Eq. ([4. 1|) is forward-backward symme- 
tric. In the continuum region of Ktt or tttt masses, there will be additional partial waves, 
as well as long-range effects associated with bremsstrahlung from B Kn and B tttt, 
followed by internal conversion of the photon. While such effects are difficult to calculate, 
it is conceivable that they lead to larger CP-violating asymmetries in the continuum region 
of -B — > Knee and B Tiiree. 
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APPENDIX A: FORM FACTORS 



In this Appendix we list the B K* and B ^ p form factors of Mehkhov and Nikitin 
using their "Set 2" of the Isgur-Scora-Grinstein-Wise parameters [|14|. (See Table [rV| .) 
TABLE IV. The B ^ K* and B ^ p form factors of Melikhov and Nikitin g. 



Form factors 



B ^ K* 



B- 



9+{q') 
9-{q^) 



2 \ -2.61 

0.048 GeV-i I 1 " ^72 



1.61 GeV 1 



+ 



5.862 7564 ) 

.2 \ -2.85 



-0.036 GeV"^ ( 1 



7.332 y 

„2 \ -3.28 



0.0037 i ^ 



-0.28 1 



6.572 _ 

.2 \ -2.62 



0.24 1 



6.672^ 

,2 \'-2.58 



6.592 



0.036 GeV"^ I 1 - ^ 
1.10 GeV I 1 



2 \ -2.75 



6.552^ 

2 4 \ -1 



+ 



5.592 7104^ 

2 N -3.04 



-0.026 GeV^i I 1 



7.292, 

,2 \ -3.42 



0.003 GeV-2 I 1 



-0.20 1 



6.432 

-.2 \ -2.76 



0.18 1 



6.572, 

,2 \ -2.73 



6.502 



APPENDIX B: THE FUNCTIONS Ii, . . . , I9 



h = (3H-siSp 



x\\xi^\' + \xn\') + {\yLV + \yRV) 



sin^ 6'p + 



The functions /i, . . . , /g appearing in the expression for /(s;, sp, 6^/, 6p, (j)), Eq. ( p.8| ), are 
jiven by 



ik-qr{\yL\' + \yRn 



+ X\\z',\' + + 2X\k . q)Re {ylz'^ + y^p) 



cos^ Op 



(Bl) 



+ X^d^l^ + + 2X^(fc ■ g)Re {ylz'^ + 



/3 = (3hiSp[x\\xL? + IxpH - {\yL? + Iz/pH} sin^ ^p, 



Sin^^p- [(fc.g)2(|^^|2+|^^|2) 

cos2^p}>, (B2) 

(B3) 
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h = (5\siSpY/^[x''Re [ylz'^ + y^z'^) + {k ■ + \yR?) } sm 2^p, (B4) 

h = 2/?2x(s,sp) V2|x2Re - + (A; • g)Re {xly^ - x^^yn) } sin 29 p, (B5) 

/e = AI3''XsiSpRe {xlyL - x%yR) sin^ Op, (B6) 

/y = 2/52x2(siSp) V^im {ylz'^ - sin 2^p, (B7) 

Is = /?'X(s,sp)^/2|x2lm + x*j,z'p) + {k ■ q)lm {xlyL + x*j,yR)'j sin 26 p, (B8) 

h = -2P^XsiSplm {xlyL + x^pyp) sin^ Op, (B9) 
where k ■ q = (M| — s/ — sp)/2. 
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